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The  apparent  inhomogeneous  nature  of  the  equatorial  electron  density  Irregu¬ 
larities  indicate  that  a  new  approach  to  the  scintillation  theory  is  required. 
Rased  on  the  available  In  situ  electron  density  measurements  and  theoretical 
lnform.it  Ion  a  simple  t wo-d Imens Iona  1 ,  deterministic  model  of  the  equatorial 
bubble  Is  developed.  The  parabolic  equation  for  a  wave  traversing  such  a  bubble 
Is  solved  numer leal lv.  The  computed  amplitude  pattern  suggests  that  wave  diffrac 
t Ion  on  edges  formed  hv  sharp  horizontal  electron  density  gradients  m;.y  be  a 
possible  cause  of  sc  1 nt 1 1 lat Ion  at  gigahertz  frequencies.  ComoarlBon  with  the, 
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result*  obtained  fur  <<  model  of  random  Irregular  it  ies  with  the  same  power 
spectrum  shows  that  the  medium  character tred  by  sharp  electron  density 
gradients  leads  to  stronger  scintillation  in  a  wide  frequency  range. 


The  amplitude  pattern,  scintillation  index, frequency  dependence  and 
scatter  plots  computed  for  the  model  agree  well  with  observations. 
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1.  INTRODUCTION 

Recent  activity  in  the  equatorial  scintillation  studies 
was  enhanced  by  the  discovery  of  seint i 1 lat inq  gigahertz 
frequency  radio  siqnals  (Skinner  et  al.,  1971)  and  by  the 
new  and  fascinating  experimental  and  theoretical  results 
for  the  irregular  structure  of  the  equatorial  ionosphere. 
Several  papers  reviewing  equatorial  scintillation  phenomena 
and  theory  and  measurements  of  equatorial  irregularities 
were  published  recently  by  Farley  et  al.,  (1970);  Rasu  and 
Kelley,  (1977);  Rasu  and  Kel ley  (1979);  Ossakow  et  al., 
(1979)).  Interested  readers  should  refer  to  these  papers 
for  more  detailed  references.  In  Section  2  only  a  brief 
summary  of  results  pertinent  to  the  future  discussion  is 
given.  As  will  become  evident,  equatorial  irregular  structure 
can  not  be  considered  as  random,  mainly  because  of  the  pres¬ 
ence  of  sharp  gradients  of  electron  concentration.  Under 
such  circumstances,  a  deterministic  approach  to  scintilla¬ 
tion  is  certainly  more  correct.  The  approach  used  in  this 
report  is  based  on  the  assumption  of  validity  of  the  para¬ 
bolic  equation  describing  the  complex  wave  amplitude  for  a 
wave  propagating  in  a  model  irregular  ionosphere.  The  model, 
which  is  deterministic,  is  based  on  the  available  satellite 
m-situ  data  and  other  exper i ment a  1  and  theoretical  informa¬ 
tion.  It  is  described  in  Section  3.  Section  4  and  the 
Appendix  discuss  the  numerical  method,  algorithm  and  computer 
program  used  to  solve  most  efficiently  the  parabolic  equation. 
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As  shown  in  Section  5,  model linq  proved  clearly  the 
importance  of  the  presence  of  sharp  horizontal  electron 
concentration  gradients  in  producing  strong  scintillation 
of  gigahertz  frequency  radio  waves.  Although  the  model  is 
deterministic,  some  statistical  parameters  of  radio  waves 
have  been  computed,  and  frequency  dependence  of  the  scintil¬ 
lation  index  was  found  to  bo  in  good  agreement  with  observa- 
t ions . 
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2.  SUMMARY  OF  CURRKNT  EXPERIMENTAL  AND  THEORETICAL 

UNDERSTANDING  OF  EQUATORIAL  IRREGULARITIES 

Most  informat ion  about  the  irregular  structure  of  the 
equatorial  ionosphere  has  come  from  back-scatter  radar  and 
in-situ  satellite  and  rocket  measurements.  Several  campaigns 
of  simultaneous  measurements  of  different  ty|H?s  provided 
deeper  insight  into  the  nature  of  equatorial  irregularities. 
Jicamarca  back-scatter  radar  data  show  that  3-m  size  irregu¬ 
larities  initially,  right  after  sunset,  occupy  relatively 
thin  layers  below  the  F-layer  peak  which  rapidly,  sometimes 
in  a  matter  of  several  minutes,  thickens  to  several  hundred 
kilometers  extending  well  above  the  F-peak  (Farley  et  al . , 
1970;  Woodman  and  La  Hoz ,  1976).  The  extended  irregular 
structures  are  called  "plumes,"  which  usually  move  with 
the  velocity  of  the  order  of  100  ms-',  occasionally  ex¬ 
ceeding  300  ms-'.  Simultaneous  VHF  scintillation  measure¬ 
ments  (Basu  ct  al.,  1977)  show  that  as  long  as  the  irregular 
layer  is  thin  scintillation  intensity  is  low  or  moderate 
and  the  fading  rate  is  small.  However,  the  stage  of  fully 
developed  plumes  is  accompanied  by  intense  (>20dB  on  VHF) 
and  rapid  scintillation.  At  times  of  intense  VHF  scintilla¬ 
tion,  qiqahertz  frequency  radio  signals  also  scintillate. 

This  scintillation  is  unexjicctedly  strong,  exceeding  8  dB 
at  1.5  GHz  and  4  dB  at  4  GHz.  Scintillation  indices  as 
large  as  0.74  at  1.5  GHz  and  0.63  at  4  GHz  were  reported 
(Taur,  1976).  Evidently,  in  such  cases,  the  weak  scattering 
scintillation  theory  is  not  valid.  Discussion  of  possible 


4 


causes  of  strong  gigahertz  frequency  scintillation  is  one 
of  the  objects  of  this  report. 

It  has  been  pointed  out  (Oasu  et  al.,  1978)  that  3-m 
and  km-size  irregularities  responsible  for  VHF  scintillation 
coexisted  only  in  the  developing  stage  of  plume,  and  in  the 
later  evening  hours  scintillation  did  not  correlate  with 
plumes.  This  led  Hasu  et  al.  to  speculate  about  the  change 
of  the  cut-off  scale  size  of  the  irregularity  spectrum. 

In-si tu  rocket  and  satellite  measurements  provided 
further  information  about  equatorial  irregularities. 

OGO-6  data  showed  that  irregularities  are  usually  asso¬ 
ciated  with  large  depletions  of  the  electron  concentration 
known  as  bubbles  (McClure  et  al.,  1977).  These  bubbles 
move  upward  and  westward  with  the  velocity  of  the  order  of 
150  ms~*  with  respect  to  the  background  plasma  (which  at 
night  rotates  eastward) .  The  observed  east-west  size  of 
the  bubbles  ranges  from  10  to  several  hundred  kilometers  and 
the  depletions  in  bubbles  up  to  three  orders  of  magnitude 
were  observed.  Figure  1  illustrates  the  irregular  struc¬ 
ture  inside  the  bubble.  The  very  sharp  gradients  and  the 
peaks  are  the  most  striking  features.  Similar  structures 
have  been  observed  below  the  F-oeak  durinq  the  rocket  flights 
(Costa  and  Kelley,  1978).  Costa  and  Kelley  (1978)  have  shown 
that  these  kinds  of  structures  have  a  one-dimensional  power 
sfiectrum  of  the  form  k-2,  resembling  the  spectrum  of  turbu¬ 
lent-like  irregularities.  From  the  point  of  view  of  the 
scint i 1 lation  theory,  this  is  a  very  important  result.  The 
irregularity  power  spectrum  plays  a  very  important  role  in 


the  statistical  theory  of  scintillation  which  makes  use  of 
the  statistical  properties  of  a  medium  to  derive  statistics 
of  a  scintillating  wave.  It  is  assumed  that  the  irregular¬ 
ities  are  random  and  the  background  medium  is  homogeneous, 
or  at  least  locally  homogeneous.  However,  in  the  presence 
of  sharp  electron  concentration  gradients  none  of  these  two 
assumptions  are  exactly  valid,  and  the  scintillation  theory 
based  on  other  models  may  be  more  appropriate  for  the  situa¬ 
tion.  Furthermore,  mechanisms  other  than  scattering  and 
diffraction  have  been  suggested  to  be  efficient  in  producing 
scintillation  (Crain  et  al.,  1979).  This  (>oint  is  discussed 
in  more  detail  later  on. 

Recent  simultaneous  total  electron  content,  airglow, 
and  scintillation  measurements  provided  further  evidence 
that  irregularities  responsible  for  scintillation  are  asso¬ 
ciated  with  depletions  in  electron  concent  rat  ion  (Koster, 

1976;  Weber  et  al.,  1978;  Yeh  et  al.,  1979).  This  asso¬ 
ciation  is  especially  clear  in  the  premidnight  period. 

Theoreticians  were  very  successful  in  explaining  the 
observed  irregular  structure  of  the  equatorial  ionosphere. 

The  hierarchy  of  plasma  instabilities  suggested  by  Haerendel 
(1974)  could  produce  a  wide  spectrum  of  observed  irregu¬ 
larities.  First  in  the  hierarchy  is  the  collisional  Rayleigh- 
Tayior  instabi 1 ity,  which  acts  on  the  steep  bottomside  gradient 
of  electron  concentration  and  produces  large-scale  verti¬ 
cally  oriented  irregularities  with  sharp  horizontal  gradients. 
Favorable  conditions  for  the  cross-field  instability  are 
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created.  This  kind  of  instability  leads  to  i rreqularities 
with  wave-vectors  oriented  vertically  and  sizes  roughly 
between  0.5  to  5  km.  When  the  amplitude  of  these  irregu¬ 
larities  is  large  enough,  the  collisionless  Ray le igh-Tay lor 
instability  becomes  active  and  produces  irregularities  with 
sizes  50  to  500  m.  Finally,  kinetic  drift  waves  can  grow 
upon  these  irregularities  after  t hoy  reach  large  amplitudes. 
Irregularities  generated  by  this  mechanism  will  have  sizes 
comparable  to  the  ion  gyro-radius  (several  meters). 

Numerical  simulation  of  the  collisional  Rayleigh- 
Taylor  instability  (c.f.  Scannapieco  and  Ossakow,  1976; 
Ossakow  et  al.,  1979)  showed  that  originally  bottomside 
plasma  depletions  can  rise  beyond  t ho  F  peak.  Figure  2 
from  Ossakow  et  al.  (1979)  shows  the  evolution  of  a  bubble. 
Percentage  depletion  as  large  as  85*  was  reached  and  the 
rise  velocity  of  t  lie  bubble  was  160  ms*  '  .  Note  the  sharp 
electron  concentration  gradient  on  the  leading  edge  of  the 
bubble.  An  analytical  model  of  uprising  bubbles  has  been 
developed  by  Ot t  (1978)  in  J he  col  1 i s i nn-domi nated  and  the 
col  1 lsionless  cases.  Kelley  and  Ott  (1978)  showed  that 
rising  bubbles  generate  a  wake  with  vortices.  Oriqinally, 
such  vortices  have  a  characteristic  size  of  the  order  of 
the  bubble  size.  In  the  presence  of  the  background  plasma 
density,  gradient  vort ices  will  mix  regions  of  differing 
density  and  create  density  irregularities  with  wavelengths 
both  larger  and  smaller  than  the  original  size  of  the  bubble. 
The  fmwor  spectrum  of  density  i rregul ar i t iet  will  be 
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proportional  to  k"*  for  irregularities  with  sues  between 
the  bubble  size  and  the  dissipation  scale  (  28  m  for  the 
case  considered  by  Kelley  and  Ott) •  This  seems  to  contradict 
the  observed  k"*  spectrum  but,  as  pointed  out  by  the  authors, 
steepening  of  the  spectrum  is  most  likely  due  to  the  sharp 
edqes  of  bubbles. 
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1.  BASIS  OF  THE  SCINT I LI.AT ION  MODEL 

The  scintillation  theory  relates  statistics  of  the  wave 
parameters  to  statistics  of  the  dielectric  permittivity 
fluctuations  «,  and,  therefore,  of  electron  concentration 
fluctuations.  If  dielectric  pormi t t i vi ty  fluctuations 
*  ,(,,2)  can  be  considered  as  a  Gaussian  random  field,  then 
statistics  ot  ,,  and  therefore  of  the  complex  amplitude, 
are  completely  determined  by  the  first  two  moments  of  the 
field:  its  mean  and  autocorrelation  function  or  structure 

function.  If  it  is  further  assumed  that  the  field  is  homo- 
aeneous,  instead  ot  the  autocorrelation  or  structure  function, 
the  |Kwr  spectrum  can  be  used.  This  assumption  is  always 
made  in  the  scintillation  theory.  It  prevents  us  from  con¬ 
sider  mu  the  fields  of  •  ,  with  sharp  qradients,  because  in 
this  ease,  •  (  is  certainly  1 nhomoqeneous .  On  the  other  hand, 
sharp  elect  ron  concent  rat  ion  qradients  are  the  most  pro¬ 
nounced  features  of  the  irroqular  equatorial  ionosphere. 
Consequently,  some  new  approach  to  the  scintillation  is 
needed.  In  addition  to  that,  sharp  discont inuit res  in  c, 
can  produce  qua  1 1 1 a t i ve 1 y  new  effects  leading  to  scintilla¬ 
tion  (Crain  et  al.,  1979). 

In  the  rest  of  this  section,  the  basis  for  the  scintilla¬ 
tion  model  free  from  assumptions  on  statistics  of  c,  is 
deve loped . 

1.1  Parabolic  Equation 

A  starting  point  for  our  consideration  is  the  so-called 
parabolic  equation  (cf.  Tatarski,  1971).  For  a  wave  propagating 
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in  the  z  direction,  it  takes  the  form 

„  Ju  , 

-2 i k  —  ♦  V_ ' u ♦ k ' i | u  -  0  (3.1) 

:>z  ' 

where  u  is  the  complex  wave  amplitude,  •  ,«»-••>  is  the  di¬ 
electric  permittivity  variation  relative  to  some  reference 
value  <t>,  V  is  the  transverse  to  the  direction  of  propa- 
qation  part  of  the  Laplacian  operator,  and  k=2n/l  is  the 
wave  number  in  the  reference  medium. 

To  apply  (3.1)  it  is  required  (Tatarski,  1971)  that: 

(A)  the  characteristic  scale  of  (  variat ions  is  much 
larqer  than  the  wavelenqth. 

( B )  wave  attenuation  due  to  t)u*  presence  of  irroqu- 
lanties  is  small  over  a  distance  equal  to  the 
wavelenqth. 

In  addition,  two  sufficient  conditions  must  be  satisfied: 

(C)  the  {>ower  reflected  back  is  neqliqible, 

(D)  the  Fresnel  approximat ion  is  valid. 

If  conditions  (A)  -  (D)  are  satisfied,  equation  (3.1) 
can  be  used  in  any  medium,  stochastic  or  deterministic. 

In  the  stochastic  medium,  the  reference  dielectric  per¬ 
mittivity  is  simply  the  averaqe  value,  while  in  the  deter¬ 
ministic  medium  it  can  be  any  value.  In  our  calculations 
*c>  has  been  chosen  to  be  equal  to  l,  which  is  the*  value 
c  assumes  outside  the  irreqular  rcqion. 

Condition  (A)  requires  iei  (  (V(  j'1  to  be  much  larqer 
than  the  wavelenqth.  Condition  (0)  is  satisfied  whenever 


(cf.  Tatarski,  1971) 
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j^TjjX  <  <  1 .  (3.2) 

and,  if  only  kt<<i,  this  inequality  is  valid  even  in  the 
far  zone  where  a  distance  L  traversed  by  a  wave  is  larger 
than  (k  i  )  . 

Conditions  (R)  and  (C)  are  much  more  difficult  to 
discuss.  To  have  some  impression  about  the  order  of  magni¬ 
tude  of  wave  attenuation  in  a  deterministic  medium,  con¬ 
sider  the  case  when  •  ,  is  a  Gaussian  function  of  x  and  y 

t  ,  (x,y)  -  exp(-p'/r?)  p-(x,y)  (3.3) 


and  this  i  rrt*gu  1  ar  l  ty  extends  throughout  a  distance  L  in 
the  direction  of  propagat ion .  Most  of  the  attenuation  re- 
sults  from  the  single  scattering,  therefore,  if  the  power 
extracted  from  the  incident  wave  by  single  scattering  is 
small  over  a  distance  equal  to  the  wavelength,  the  condi¬ 
tion  (B)  is  satisfied.  Assuming  that  the  Fresnel  approxi¬ 
mation  is  valid,  a  single  scattered  wave  field  is  given  by 
the  equation: 


u,  ( .> ,  L ) 


*  u»  (  dz*  ff 

4-  1  L-»S*  LI 

0 


(p*  ,  z'  lexpl- 


*k<P-P’)  , 
2 (L-z  * ) 


dp ' 


(3.4) 


where  u,  is  constant  for  a  plane  incident  wave.  For  sim¬ 
plicity,  we  set  o*0,  therefore,  the  up|>er  estimate  of  the 
scattered  power  is  obtained.  Substituting  (3.3)  into  (3.4), 
integrating  and  taking  the  absolute  value  of  u,  leads  to 
the  result 
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l  ^  r  ^  2 

k  _L ' « 
16 


(tn7p«V)  . 


1 1 


(3.5) 


P*’  =  1  ♦  (2L/r’k)',  $  «  arctan  (2L/r‘  k)  , 

which  is  the  fraction  of  the  incident  jwwor  attributed  to 
a  sinqle  scattering.  The  inverse  of  the  extinction  length 
d  is  obtained  by  dividing  (3.5)  by  t he  slab  thickness  L. 
Extinction  length  is  a  measure  of  the  wave  attenuation. 

If  d  */  k- •  1  condition  (B)  is  satisfied. 

Consider  first  the  case  2J./kr  ’  1  .  It  is  easy  to  show 

that  in  this  case 


d”l/k 


kb.  ,  •  /4  . 
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(3.6) 


If  (3.6)  is  multiplied  anti  divided  by  2kr',  the  condition 
(B)  can  be  written  in  t he  following  form 


k '  r  ‘ 


i  o 


21. 

kr‘ 


<<  1  . 


(3.7) 


Note  that  an  irregularity  of  size  r  causes  the  phase 
change  kr*  »#/2.  Therefore,  inequality  (3.7)  means  that  for 
2L/kr*<<l,  the  condition  (B)  is  satisfied  even  for  large 
phase  variations  introduced  by  a  single  irregularity. 

For  2L/kr'>>l  the  value  of  d~*/k  can  bo  approximated 
by  the  following  expression: 

ii  i  fn(2L/kr‘) 

d  */k  -  -■  ---•  |fn(2I./kr?)  )2  =»  kJL?i,#*  | — — — ; - 1 


2L/kr? 


(3.8) 


and  the  condition  d-1/k-<l  means  that  the  total  phase  chanqe 
kL*  , #  can  not  be  exceedingly  large. 
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To  discuss  condition  (C)  again,  assume  the  Gaussian 
irregularity.  One  can  expect  that,  whenever  the  scattering 
angle  is  small  not  too  much  power  is  reflected  back.  To 
estimate  the  scattering  angle  note  that  in  the  plane 
y=const  it  can  be  expressed  in  terms  of  the  phase  $  by 

.-lit  (3.9 

k  .lx 

where  $-k  ,L/2.  In  our  irregularity  model,  the  scattering 
angle  is  maximum  for  y=0  and  x^r/*'?,  where  it  assumes  the 
va  1  up 


I«i  |  rt  kL.  |  q 

«  v  0.46  1  !  0.46 - 11 

ma  x  ,  k  r 


(3.10) 


Condition  (C)  is  satisfied  provided  ‘max<r^-  It  is  clearly 
seen  from  (1.10)  that  with  kr->l  this  puts  a  limit  on  the 
total  phase  change. 

With  strong  irregularities  filling  up  the  extended 
bubble  structures,  assumptions  leading  to  the  parabolic 
equation  could  be  violated  especially  at  lower  frequencies. 

In  fact,  condition  (R)  is  violated  in  our  model  calculations 
when  scintillation  at  a  frequency  of  116  MHz  is  considered 
to  l>e  due  to  irregularities  at  the  center  of  the  well- 
developed  bubble  (Section  3.2).  A  rough  estimate  shows 
that  in  this  case  d“*/k  is  larger  than  one  for  irregularities 
with  the  scale  size  sm.il  ler  than  200  m. 


1.2  Irregularity  model 

A  parabol ic  equat ion  can  be  solved  numerically  only  if 
•  ,(e,z)  is  known  at  each  point  of  a  three-dimensional  space. 


1  3 

Because  of  the  direct  proportionality  to  electron  concen¬ 
tration  fluctuations  AN,  i  ,  can  be  determined  provided  AN 
is  known  from  measurements.  Unfortunately,  satellite  and 
rocket  in-situ  measurements  give  information  only  about  AN 
along  the  satellite  (rocket)  path.  Any  attempt  to  obtain 
from  these  data  a  three-dimensional  structure  of  irregu¬ 
larities  would  involve  assumptions  about  their  statistics. 

This  is  feasible  if  the  data  show  randomness  of  AN.  Figure  1 
clearly  shows  that  this  is  not  the  case:  in  addition  to 
the  large  scale  trend,  sharp  gradients  of  AN  are  present 
inside  the  bubble. 

The  two-dimensional  model  adopted  lor  our  calculation 
is  very  simple.  It  is  assumed  that  at  each  height  the  hori¬ 
zontal  relative  variations  of  electron  concentration  AN/N 
are  the  same  as  those  displayed  in  Figure  1.  The  concen¬ 
tration  at  the  bubble's  edge  is  taken  as  a  reference.  Allow¬ 
ance  is  made  for  its  variations  with  the  height  and,  in 
particular,  the  parabolic  profile  is  adopted. 

N  ( h )  *  N#  (1- (h-h#)  VH?  I  (3.2.1) 

where  N,  is  the  maximum  electron  concentration  at  the  heiqht 
h,  and  H  is  the  scale  height.  In  most  cases  we  set  h„=350  km, 
No  =  2.S*10C’  electrons  cm-'  and  H  110  km. 

As  was  noticed  in  Section  2  theoretical  simulations 
show  that  in  the  course  of  the  bubble  development  its 
leading  edge  becomes  sharper  while  the  trailing  edge  seems 


to  be  more  blunt.  To  take  this  into  account  (3.2.1)  is 
multiplied  by  some  weighting  function.  The  weiqhting 
function  corresponding  to  what  we  call  initial  stage  is 
shown  in  Figure  3a,  while  Fiqure  3b  displays  the  weighting 
function  for  the  developed  staqe  of  the  bubble. 

The  model,  is  summarized  in  Figures  4a  and  4b,  which 
show  electron  concentration  variations  AN  at  several  levels 
inside  the  bubble  for  the  initial  and  developed  stages, 
respectively. 

Our  model  is  far  from  being  perfect.  In  particular 
it  does  not  take  into  account  small-scale  vertical  varia¬ 
tions  of  the  electron  concentration.  Hut  at  least  it  is 


free  from  any  assumptions  about  the  statistics  of  AN  and 
should  correctly  answer  the  question  about  the  role  of 
vertically  extended  sh.jrp  gradients  in  producing  scintillation 


4. 
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NUMERICAL  METHOD  FOR  SOLVING  PARABOLIC  EQUATION 
The  parabolic  equation  (3.1)  can  bo  solved  analyti¬ 


cally  only  for  a  very  tew  simple  forms  of  the  function 
Ci(p,z).  A  general  numerical  approach  is  needed.  We  are 
interested  in  solving  (3.1)  when  the  dielectric  permittivity 
is  a  function  of  x  and  z  only,  thus,  V_,=3,/3xJ.  For  a 
plane  incident  wave  the  initial  value  for  this  equation 
is  u  ( x ,  0 ) -const*  1  .  The  change  of  variables  t  kz,  f.  =  kx  leads 
to  the  equation 

*u  i  a*'u  i 

♦  ♦  .,(*.,  t )  u  0  (4.1) 

>t  2  *  2 

Basically,  two  methods  for  the  numerical  treatment  of 
parabolic  equations  exist:  explicit  and  implicit  method 
(cf.  Ames,  1969).  Although  t  hi1  explicit  method  is  very 

simple  and  does  not  involve  the  troublesome  solution  of  a 

large  set  of  equations,  it  gives  a  somewhat  imperfect  model 
for  a  parabolic  equation  ( Arm's ,  1969),  and  places  a  stringent 
restriction  on  the  step-size  in  the  t-direction.  Because  of 
this,  we  have  chosen  as  an  implicit  scheme  a  special  case 
which  is  the  well-known  Crank-N ichol son  difference  scheme. 

A  thorough  discussion  of  the  stability  condition,  trunca¬ 
tion  error  and  convergence  has  been  made  by  Crandall  (1955) 

(see  also  Ames,  1969)  for  the  diffusion  equation  and  for  a 
special  case  of  the  equation  analogous  to  (4.1)  by  Rorodz iewicz 
(1978).  The  scheme  makes  use  of  the  following  analog  to  the 
second  derivative: 


u 
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no  corresponding  relation  was  derived.  However,  computer 
runs  showed  that  a  satisfactory  compromise  between  the 
error  and  convergence  is  attained  if  o=0.51  and  0.5<r<2.5. 

Substitution  of  (4.2),  (4.3)  and  (4.4)  into  the  para¬ 

bolic  equation  (4.1)  leads  to  the  following  difference 
equat ion : 


l  T 
2h‘ 


ou 


n* 1 
k  ♦  1 


i  r 

it  n 

n  ♦  1  i  l 

n*  1 

+  (1- 

-r  0  ♦ 

—  of.  ) 

ou  .  = 

h 

2  k 

k  2h^ 

k-  1 

it 

n 

i  t 

it  n 

n 

"  h* 

( l-o  )  u.  ,-*-11  ♦ 
k*  1 

h”  <1-’)  - 

2-  <l-°>fk! 

|uk 

it  n 

- ; ( 1 -o )  (4.6) 

2h '  k-  1  ’ 


k  »  ♦ 1 ,  *2,  ...  ♦  K 

n  *  0,  1,  2,...  N 

n  n  ♦  1  n 

whore  fk-(c,k 

In  the  discrete  scheme  the  initial  value  u(x,0)=l  is 
written  u^*l.  Because  of  the  assumed  symmetry  of  t  ,(x,z) 
with  respect  to  the  z  axis,  un  and  the  system  (4.6)  of 

2K  equations  reduced  to  the  system  of  K  equations.  It  con¬ 
tains  K*2  unknowns  un  with  K  unknowns  of  interest  to  us 

k 

(k«l,2,  ...  ,  K) ,  and  two  unknowns  u™  and  u^  ^  beyond  the 
boundaries  of  the  considered  region.  The  usual  procedure 
is  to  introduce  the  so-called  "false"  boundaries  at  k*0 
and  k»K*l.  If  the  derivatives  of  u  at  the  "true"  boundaries 


are  known  these  "false"  boundaries  can  be  easily  derived. 


18 


Wo  sot  un=0  and  u*'  =0  which  introduces  some  discontinuities. 

0  k  +  1 

As  noticed  by  Ames  (1969),  if  the  stable  scheme  is  used,  the 
effect  of  such  discontinuities  does  not  penetrate  deeply 
into  the  field  of  integration.  In  this  way,  the  system  (4.6) 
becomes  tridiagonal  with  K  unknowns.  It  can  be  written  in 
the  matrix  form 

BnUn+l  -  (Bn-tAn)Un,  n=0 , 1 , 2 . N  (4.7) 


whore  is  the  column  matrix  with  u|^  as  elements.  Matrices 
Am  and  Bn  are  tridiagonal  with  the  following  elements: 


n 

A 

k  ,  k 


n  n 

A  -  A 

k  ♦  1  ,  k  k  ,  k  ♦  1 


n 

A  0  for  k/  i  k/  ]  ♦ 1 
k  ) 


k  =  1 , 2 . K 

(4.8) 


n  1  *  L-  n  n 

B  -  1  - 1  i  n  (  -  —  )  B  B 

k.k  h‘  2  k  ♦  1  ,  k  k,kU 


1  i  o 

2h^ 


k  -  1  ,  2 


B"  -0  for  k/i  k/)*l 
k  *  i 


,K 

(4.9) 


The  solution  of  the  matrix  equation  (4.7)  is 
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n+  1 


n  n  ~  *  n  n 
U  -  (B  )  i  (A  U  ) 


(4.10) 


To  get  the  second  term  on  the  right-hand  side  of  (4.10),  An 
is  first  multiplied  by  Un ,  already  known  from  the  previous 
step,  and  then  the  matrix  equation  BnX~AnUn  is  solved. 
Because  Bn  is  tridiagonal,  the  matrix  inversion  can  be 


avoided  by  usinq  the  Gaussian  elimination.  The  so-called 
Thomas  algorithm  is  most  often  used.  The  system  of  equa¬ 
tions  is  written  in  the  form 


Notinq  that  x  the  unknowns  xu  are  successively 

found  from  (4.12).  It  has  been  shown  that  this  scheme  does 
not  introduce  larqe  errors  due  to  round-off  errors  in  the 
calculations  (Ames,  1969)  and  is  very  efficient  in  terms 
of  computer  time. 

A  complete  proqram  for  solving  the  parabolic  equation 
is  given  in  the  Appendix. 


S.  NUMERICAL  RESULTS  ANl)  DISCUSSION 

The  parabol ic  equation  (3.1)  has  been  solved  for  a 
model  equatorial  plasma  bubble  filled  with  1 r reqular i t ies 
whose  horizontal  structure  is  the  same  as  the  one  measured 
in-situ  by  McClure  et  al.  (1977)  and  shown  in  their  Figure  7. 
The  satellite  altitude  was  approximately  280  km  and  probably 
did  not  change  appreciably  during  the  passage  through  a 
bubble  whose  horizontal  size  was  ^50  km.  The  F-peak  height 
at  the  time  of  the  measurements  was  approximately  400  km. 

After  enlargement,  the  curve  in  Figure  7  of  McClure 
et  al.  was  digitized  every  42  m  which  is  approximately  10  m 
less  than  the  sampling  interval  of  in-situ  measurements . 

This  sets  the  lower  limit  on  t he  irregularity  scale-size 
or  the  scale  length  of  the  electron  density  gradient.  It 
is  still  much  less  than  the  Fresnel  zone  dimension  for  a 
typical  distance  between  the  bottom-side  of  the  irregular 
layer  and  the  receiver  (200-100  km)  and  the  highest  wave 
frequency  considered  in  this  report ,  which  is  6  GHz. 

To  have  some  impression  about  the  irregularity  spec¬ 
trum,  the  digitized  data  were  high-pass  filtered  to  remove 
any  components  larger  than  10  km  and  the  power  spectrum 
was  computed.  Figure  S  shows  that  the  spectral  ampli¬ 
tude  roughly  follows  the  ►  “1  law  except  for  the  extreme 
point  at  the  low  frequency  part  of  the  spectrum.  This  is 
a  very  typical  spcct rum  of  ionospheric  F  region  irregularities 
(Dyson  et  al.,  1974). 
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We  assumed  that  the  rel.it  1  ve  electron  density  fluctua¬ 
tions  at  each  height  inside  a  bubble  follow  the  same  hori¬ 
zontal  pattern  as  one  given  in  Figure  1.  The  background 
density  is  parabolic  as  described  in  Section  2.  With  the 
weighting  function  corresponding  to  the  "initial  stage" 
of  a  bubble,  before  it  rose  above  the  F-peak,  the  horizontal 
absolute  electron  density  variation  AN  at  four  different 
levels  inside  the  bubble  is  given  in  Figure  4a.  The  largest 
AN  *1.7*10’  electrons • cm” 5  is  reached  close  to  the  center  of 
the  bubble.  If  the  bubble  thickness  L  is  defined  as  a  distance 
between  levels  at  which  the  weighting  function  falls  to  e 
we  have  L“30  km. 

The  second  extreme  case  to  be  discussed  corresponds 
to  the  "developed  stage"  of  a  bubble  (Figure  4b).  In  this 

C, . 

case,  the  maximum  ’N  is  approx i ma te 1 y  egual  to  2*10  elec- 
trons*cm  '  and  the  bubble  thickness  I.  1  1  r>  km. 

The  scintillation  pattern  resulting  from  our  com¬ 
putations  for  a  wave  of  frequency  1.5  GHz  passing  through 
th?  initial  stage  bubble  is  shown  in  Figure  k.  The  real 
and  imaginary  comjionents  of  the  complex  amplitude  as  well 
as  amplitude  anti  phase  in  r  units  are  depicted.  The  most 
interesting  feature  is  the  correspondence  between  the  position 
of  strong  amplitude  outbursts  and  the  location  of  sharp  elec¬ 
tron  density  gradients  in  Figure  1.  This  seems  to  suggest 
the  refractive  scattering  mechanism  proposed  by  Crain  et  al. 
(1979)  as  a  cause  of  appreciable  scintillation.  However, 
as  we  show  later  in  this  section  the  diffraction  on  edges 
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more  easily  explains  the  obtained  feature.  It  is  also 
interest  inq  to  note  one  more  striking  p<-  liarity,  namely, 
the  weak  dependence  of  the  amplitude  outbursts  on  the 
absolute  value  of  electron  density  deviation  AN. 

For  a  develop'd,  thick  bubble  the  edqe  effect  is  largely 
obscured  by  amplitude  fluctuations  attributed  to  scattering 
or  diffraction  by  smoother  structure  (Figure  7).  Neverthe¬ 
less,  stronger  amplitude  scint i 1 lat ion  can  still  be  noted 
close  to  sharp  gradients  of  AN.  On  lower  frequencies  the 
edge  effect  can  hardly  be  distinguished  from  scattering 
(Figure  8b).  However,  its  presence  can  be  seen  from 
Figure  8a,  which  shows  the  amplitude  pattern  at  300  km 
height,  i.e.,  inside  the  bubble.  Fvidently,  the  edge 
diffraction  develops  amplitude  fluctuations  at  shorter 
distances  than  does  the  scattering. 

To  discuss  further  the  effect  of  an  isolated  sharp 
gradient  in  electron  density,  we  assumed  a  wave  passing 
through  a  one-dimensional  screen  which  at  the  point  x=x0 

changes  the  wave’s  phase  from  0  to  some  constant  value  4  . 

0 

The  complex  amplitude  at  a  distance  z  is  given  by 

_  <• 

U(x,7. )  *  /  *  [  u  (  x  ’  ,  0)  exp  |  ik  i*Z*JLL\dx\ 

2  *  i  z  I  2  7. 


where 


u(x,0)  *  ol*»  for  x*x# 
*  1  for  X'X# 
Simple  integration  and  separation  of 


(5.2) 

the  real  and  imagi¬ 


nary  parts  of  u(x,z)  lead  to  the  following  formulae 


Re  u ( x , z )  *  cos :  ’ 0 

2 


2  » 


-  !  (  11-/3  cos  ( J  -  ♦,)  |C((1) 

♦  11-/2  sin  (~  -  i#))S(i)l,  (5.3) 

Im  u(x,z)  *  1  sin  A  -  -  (|l-/2  cos  (  -  4  )|S(u) 

2  0  2  4  9 

-  11-/3  sin  (J  -  ♦.  )  1 C  (.«)),  ( S .  4  ) 

where  S(i)  and  C(»)  are  the  Fresnel  integrals  of  the  argument 
i= lx-Xj ) /fc/2z .  Kquatians  (5.3)  and  (5.4)  were  used  to  obtain 
the  diffraction  pattern  at  a  distance  z  150  km  from  a  screen 
for  frequencies  of  ifeO  MHz  and  4  GHz  with  4  1  <» .  4  and  1.48, 
respectively.  The  sudden  phase  jump  was  put  at  x  -  3  km.  Re¬ 
sults  of  the  calculations  are  shown  in  Fiqure  .  A  typical 
oscillatory  pattern  can  be  noticed.  Oscillations  are  much 
faster  at  4  GHz  and  dumped  at  shorter  distances  I rom  the  dis¬ 
continuity,  therefore,  the  edqe  effect  is  more  localized  in 
this  case  than  at  lower  frequency.  Tins  explains  why  at 
hiqher  frequencies  it  is  easier  to  find  the  correspondence 
between  sharp  electron  density  qradients  and  the  character¬ 
istic  wave  amplitude  pattern.  At  lower  frequencies  the 
effect  of  closely  spaced  sharp  qradients  leads  to  a  rather 
complicated  picture  not  resemblinq  the  pattern  produced  by 
localized  discontinuities.  Thus,  the  edqe  effect  is  not 
seen  in  Figure  8b. 

Crain  et  al.  (1979)  suqqested  the  refractive  scattering 
mechanism  as  a  possible  cause  of  scintillation  at  qiqahertz 
frequencies.  This  mechanism  involves  a  total  internal 


reflection  from  vertically  extended  horizontal  discontinuities 
of  electron  density.  Under  some  strict  conditions  imposed 
on  the  slope  of  discontinuity,  the  direct  wave  interferes 
with  the  unattenuated,  reflected  wave  leading  to  stronq 
amplitude  f luctuat ions .  A  thorough  discussion  of  the  mechanism 
would  require  ray  tracing  in  the  model  bubble  with  sharp  boun¬ 
daries.  However,  if  the  maximum  grazing  angle  is  small  the 
parabolic  equation  offers  an  adequate  mathematical  description 
of  the  effect.  Thus,  we  might  expect  that  our  calculations 
take  into  account  refractive  scattering.  To  check  the  lm- 
portance  of  refractive  scattering  relative  to  the  diffrac¬ 
tion  on  edges,  we  computed  the  diffraction  pattern  produced 
by  a  screen  which  changes  the  phase  of  the  wave  proportionally 
to  ’N  given  by  tn-situ  data.  The  proport ional ity  constant 
w.is  chosen  such  as  to  achieve  the  best  agreement  between  the 
phase  variations  at  a  distance  of  ISO  km  from  a  screen  with 
those  at  a  distance  400  km  from  t he  top  of  a  bubble  when  the 
full  mcxlel  is  used.  We  wish  to  note  at  this  point  that  the 
phase  screen  was  placed  not  in  the  middle  of  the  scattering 
layer  (approx i mat e 1 y  100  km  height)  but  above  the  layer 
(cf.  Figure  la).  This  was  necessary  in  order  to  make  allow¬ 
ance  for  the  development  of  amplitude  sci  nt  i  1  lat.  i  on  as  strong 
as  that  given  by  the  solution  of  the  parabolic  equation. 

This  remark  might  be  of  some  importance  to  those  who  use 
the  strong  phase  screen  theory  in  interpreting  the  data. 

Figures  10  and  11  present  the  comparison  of  the  ampli¬ 
tude  and  phase  variations  at  160  and  4000  MHz.  as  produced 
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by  the  phase  screen  and  the  initial  stage  model  bubble. 

Close  inspection  of  the  figures  indicates  a  remarkable 
similiarity  between  amplitude  anti  phase  structures  with 
only  slightly  deeper  fading  at  )(>0  MHz  produced  by  the 
bubble.  We  conclude  that  at  least  to  the  order  of  approxi¬ 
mation  given  by  t  hi*  parabolic  equation  the  refractive 
scattering  is  either  of  no  importance  as  compared  to  the 
edqe  diffraction  effect  or  both  effects  are  indistinguishable. 

In  Figures  12  and  13  we  present  the  amplitude  fluctua¬ 
tions  at  several  frequencies  for  the  initial  and  developed 
stages  of  the  bubble,  respectively.  The  edge  effect  is 
clearly  seen  on  gigahertz  frequencies  in  the  developed 
stage  and  even  at  800  MHz  in  the  case  of  the  initial  bubble 
stage.  It  is  interesting  to  note  an  increase  in  the  fading 
rate  as  the  bubble  thickens.  This  is  clearly  seen  at  all 
frequencies  except  at  4  GHz  which  is  dominated  by  the  edge 
effect.  This  agrees  very  well  with  measurements  (Basu  et 
al.,  1977)  showing  rising  amplitude  fading  rate  when  the 
bubble  develops  and  extends  above  the  F-peak.  At  the  out¬ 
skirts  of  the  bubble  (distance  larger  than  15  km),  scintilla¬ 
tion  intensity  and  fading  rate  are  smaller  than  those  closer 
to  the  center,  mainly  due  to  the  deficiency  of  small  scale 
electron  density  irregularities.  The  fading  rate,  which  is 
also  changing  with  frequency,  is  usually  larger  at  lower 
frequencies.  This  is  in  agreement  with  observations  by 
Rino  et  al.,  (197R)  but  contrary  to  the  weak  scatter 
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scintillation  theory  which  predicts  the  dominant  fading 
period  proportional  to  the  Fresnel  zone  size  /Tz.  However, 
when  scintillation  is  stronq,  as  in  our  case  at  frequencies 
below  4  C5Hz,  the  multiple  scattering  effects  reduce  the 
characteristic  size  of  the  wave  field  fluctuations  leading 
to  the  increase  of  the  fading  rate  (Yeh  et  al.,  1975).  The 
practical  importance  of  the  scintillation  index  justifies 
its  calculations  in  spite  of  the  obvious  nonstat ionar ity  of 
the  computed  amplitude  pattern  noticeable  especially  at 
higher  frequencies.  The  S .  scintillation  index  was  calcu¬ 
lated  only  for  the  segment  of  data  from  x  0  to  15  km.  The 
segment  of  x= 1 5  to  25  km  shows  clearly  a  different  amplitude 
structure  due  to  the  lack  of  small-scale  electron  density 
l r regular i t ies .  Between  x  0  and  15  km,  scintillation  indices 
calculated  for  each  5  km  data  segment  differed  not  more  than 
25*  which  give*?  the  impression  about  the  stationarity  assump¬ 
tion.  The  frequency  dependence  of  sc int i 1 lat ion  index  is 
shown  in  Fiqure  14.  For  moderately  stronq  scintillation 
(S4  0.5)  ,  S4  is  projjort  ional  to  f”  with  v  between  1  and  1.5. 
However,  the  index  v  is  closer  to  1.5  for  the  scintillation 
index  smaller  than  0.2.  As  one  would  expect,  the  stronger 
t he  scintillation  the  weaker  is  its  frequency  dependence. 

We  must  f>oint  out,  however,  that  at  low  frequencies  the 
assumptions  leading  to  the  parabolic  equation  might  break 
down.  In  fact,  at  136  MHz  in  the  case  of  a  thick,  developed 
bubble,  the  average  intensity  is  about  15*  lower  than  the 
intensity  of  the  incident  wave.  This  means  that  15*  of  the 
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incident  power  is  backscattered  from  the  irregular  layer, 
which  no  longer  is  small.  In  other  cases  under  considera¬ 
tion,  the  backscattered  power  was  less  than  10%  of  the  inci¬ 
dent  power.  Figure  15  shows  the  case  when  the  average  signal 
level  close  to  the  center  of  the  bubble  is  much  below  the 
level  of  the  incident  wave.  In  this  particular  example, 
only  about  77%  of  the  fxiwer  is  transmitted  through  the 
bubble.  The  result  shown  in  Figure  15  was  for  frequency 
wave  of  136  MHz  passing  through  the  extended  bubble  and  the 
F-peak  electron  density  was  3.5*  10r’  electrons  •cm”  * .  It 
would  be  interesting  to  examine  the  ohservat iona 1  data  to 
find  cases  where  strong  scintillation  is  accompanied  by  a 
decrease  in  the  signal  level.  If  such  cases  could  be  found, 
this  would  indicate  the  need  for  modification  of  the  scin¬ 
tillation  theory  based  on  the  parabolic  equation. 

It  is  interesting  to  compare  the  scintillation  pattern 
produced  by  a  bubble  with  sharp  electron  density  gradients 
to  that  given  by  a  bubble  with  random  irregularities.  To 
do  that  the  electron  density  profile  shown  in  Figure  1  was 
low-pass  filtered  with  the  cut-off  wave  number  0.1  km" * . 

The  high-f roquency  comjxinent  of  the  data  was  obtained  by 
subtracting  low-frequency  part  from  the  original  data. 

The  high-frequency  congjoncnt  of  the  data  was  then  Fourier 
analyzed  and  the  amplitude  and  phase  of  each  Fourier  com¬ 
ponent  was  determined.  Following  the  procedure  used  by 
Costa  and  Kelley  (1978)  a  random  phase  angle  was  added  to 
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each  phase,  the  data  were  then  transformed  back  to  the  spatial 
domain  and  added  to  the  low-frequency  component  of  the  oriqinal 
data.  As  a  result  we  obtained  the  electron  density  profile 
with  random  h iqh- f requency  component  which  has  the  same  power 
spectrum  and  variance  of  the  electron  density  as  the  high-fre¬ 
quency  component  of  the  oriqinal  data.  Two  examples  of  such 
profile  qenerated  using  different  random  numbers  are  shown  in 
Figures  16a  and  16b.  Note  that  sharp  spikes  and  gradients  clea 
seen  in  Figure  1  are  not  present  on  these  profiles.  On  the 
other  hand,  they  contain  more  weak,  small-scale  structures. 

The  scintillation  pattern  produced  by  the  "initial 
stage"  bubble  with  horizontal  electron  density  profiles  as 
shown  in  Figures  16  was  calculated  leading  to  the  result 
shown  in  Figure  17.  At  4  GHz  almost  no  scintillation  is 
seen.  The  weakness  of  scintillation  is  confirmed  by  the 
calculated  scintillation  index  of  0.02  for  both  runs.  This 
is  approximately  half  of  the  scintillation  index  calculated 
when  the  original  profile  is  used.  The  variance  of  the 
high-frequency  electron  density  comj>onent  is  equal  to  6.6*10^ 
elect rons •cm”  1  which  is  only  slightly  smaller  than  that  for 
the  oriqinal  data  (7.7*10?  e lect rons • cm- ^ ) .  Thus  the  two¬ 
fold  decrease  in  the  scintillation  index  can  not  be  explained 
by  the  difference  in  <(AN)*>  since  in  the  weak  scattering 
scintillation  theory  the  scintillation  index  is  directly 
proportional  to  <  (AN)  Therefore,  random  electron  den¬ 

sity  fluctuations  would  have  to  be  stronqer  (in  our  case  two 
times  stronqer)  to  give  scintillation  as  strong  as  that 
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produced  by  a  medium  dominated  by  sharp  electron  density 
gradients.  Also  at  lower  wave  frequencies  the  scintilla¬ 
tion  index  is  enhanced  by  the  presence  of  sharp  electron 
density  gradients.  For  instance,  at  360  MHz  for  the  original 
electron  density  profile  S^=0.64  while  the  randomized  pro¬ 
files  lead  to  S4=0.45  and  0.52,  for  run  A  and  B,  respectively. 
Obviously,  the  difference  between  the  scintillation  indices 
is  smaller  than  at  higher  frequencies  because  of  the  effect 
of  stronq  scattering  which  reduces  the  scintillation  index 
dependence  on  <AN?>. 

Sometimes  observat iona 1  data  are  represented  in  the 
form  of  scatter  plots  on  the  plane  Imu-Reu.  For  complete¬ 
ness,  we  show  in  Figure  18  several  such  plots,  all  of  them 
obtained  using  the  same  "initial  stage"  bubble  model  and 
differing  in  the  wave  frequency  only.  No  filtering  to 
separate  the  focus  and  scatter  components  was  applied 
mainly  because  of  the  nonst a t iona r i t y  of  the  imaginary  and 
real  parts  of  the  complex  amplitude  and  difficulty  in  choosing 
the  appropriate  cut-off  wave  number.  Nevertheless,  it  is 
interesting  to  note  the  obvious  similarity  between  the 
model  plots  and  plots  derived  from  the  scintillation  data 
(cf.  Fremouw  et  al.,  1978).  Both  the  model  and  data  show  a 
clear  change  from  the  phase  dominated  scintillation  at  high 
frequencies  to  the  dominance  of  the  scatter  component  with 
the  power  evenly  distributed  between  Imu  and  Reu  at  lower 
frequencies. 

In  our  model  we  have  assumed  that  the  relative  electron 
density  fluctuation  AN/N  does  not  depend  on  the  height  inside 
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a  bubble.  This  assumption  leads  to  the  exaggerated  phase 
fluctuations  and  possibly  also  to  the  overestimated  scattering 
effects  in  amplitude.  However,  the  edge  diffraction  effect 
is  correctly  predicted. 
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6.  CONCLUSIONS 

Solution  of  the  parabolic  equation  for  the  complex 
amplitude  of  a  wave  passing  through  a  simple  model  of  an 
equatorial  bubble  revealed  the  imjKjrtance  of  sharp  hori¬ 
zontal  gradients  of  the  electron  density  in  producing  scin¬ 
tillation  effects.  Preliminary  discussion  of  the  effects 
suggests  the  diffraction  on  edqes  formed  by  sharp  gradients 
as  a  possible  mechanism  of  scintillation.  The  edge  effect 
is  particularly  important  at  higher  gigahertz  frequencies 
where  it  is  not  obscured  by  scattering.  The  important  fact 
is  that  the  edge  diffraction  leads  to  appreciable  scintilla¬ 
tion  even  if  AN  is  quite  moderate  and  L  is  small. 

The  computed  amplitude  scintillation  pattern  is  in  good 
agreement  with  the  observations,  which  indicate  an  increase 
of  the  scint i 1 lat ion  rate  with  the  increase  of  the  bubble 
thickness  and  show  a  change  in  the  scintillation  rate  from 
slow,  close  to  the  edqes  of  the  bubble,  to  fast,  at  its 
center.  The  frequency  dependence  of  scintillation  is  of  the 
type  f"  with  v  between  1  and  1.5  at  gigahertz  frequencies. 
There  is  a  tendency  toward  smaller  values  of  v  as  the  bubble 
thickens . 

The  bubble  model  we  have  used  is  deterministic.  Due 
to  the  obvious  nonhomogeneity  of  electron  density  fluctua¬ 
tions,  it  is  difficult  to  compare  the  obtained  results  with 
those  based  on  the  statistical  scintillation  theory.  Never¬ 
theless,  we  have  shown  that  the  model  of  random  irregularities 


32 


characterized  by  the  same  power  spectrum  as  the  medium  with 
sharp  electron  density  qradients  produces  less  scintillation 
not  only  at  GHz  frequencies  but  also  at  lower  frequencies 
where  the  scatterinq  effects  should  be  dominatinq.  Apparently, 
the  increased  strenqth  of  scintillation  can  be  attributed  to 
the  diffraction  on  sharp  edqes--the  scintillation  mechanism 
which  is  not  taken  into  account  by  the  statistical  scintilla¬ 
tion  theory. 
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APPENDIX 

In  this  appendix  the  algorithm  and  computer  program 
(in  CDC  FORTRAN  Extended)  used  to  solve  the  two-dimensional 
parabolic  equation  (4.1)  are  given.  The  employed  method  is 
described  in  Section  4. 

The  proqram  was  made  versatile  to  be  used  with  an 
arbitrary  model  of  the  irregular  ionosphere,  but  also  to 
make  possible  calculations  for  a  wave  which  propagates  in 
a  free  space  after  traversing  an  irregular  layer  or  a  wave 
disturbed  by  a  phase  changing  screen. 

The  input  parameters  are: 

K  -  n umbi' r  of  steps  in  the  horizontal  x  direction, 

N  -  number  of  steps  in  the  vertical  z  direction, 

TAU  -  step-size  in  the  z  direction  (in  m) , 

H  -  step-size  in  the  x  direction  (in  m) , 

Sli'iM  -  weighting  factor  (our  choise  was  SIGM*0.51), 

FRKO  -  wave  f roquency  (in  MHz), 

N 1  -  Nl’l  for  the  phase  screen  model, 

KSKir  -  number  of  steps  in  the  z  direction  recorded 

on  the  file  PASS  (which  is  the  renamed  file 
RESULT)  over  which  one  has  to  skip  to  get 
to  the  z  level  from  which  further  computations 
will  be  made.  This  is  used  whenever  one  wants 
to  continue  computations  with  the  changed  step- 
size  TAU,  or  when  further  computations  are  to 
be  made  outside  a  scattering  medium. 

XM  -  only  every  XM  step  in  the  z  direction  is  re¬ 

corded  on  files  RESULT  and  TERM 

KK  -  only  every  KK  step  in  the  x  direction  is  re¬ 

corded  on  file  TERM. 
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IC  -  used  as  an  indicator.  If  IC*1  computations 

are  made  inside  the  scattering  medium,  if  IC=0, 
for  computations  outside  a  medium  and  for  the 
phase  screen  model. 

KL  -  total  number  of  steps  in  the  i  direction  cal¬ 

culated  in  the  previous  runs.  KL  is  used  to 
calculate  the  overall  distance  passed  by  a 
wave . 

NP  -  used  as  the  reading  indicator.  If  NP^O  input 

parameters  recorded  on  the  file  PASS  are  read. 
Otherwise  it  is  assumed  that  they  were  not  re¬ 
corded  . 

K2  -  K2*0  if  the  input  file  KLFL  must  be  read.  It 

does  not  need  to  be  read  for  computations  out¬ 
side  a  scattering  medium,  but  must  be  read  at 
the  beginning  of  the  phase  screen  computations. 

PHL0  -  phase  variance  or  a  constant  converting  data 
stored  on  KLFL  to  phase  variations. 

The  above  parameters  are  universal  and  must  be  supplied 
whenever  the  program  is  used.  Below,  the  described  parameters 
which  apply  only  to  the  model  of  a  plasma  bubble  discussed  in 
the  report: 

AI.  and  7.0  -  weighting  function  at  the  top  of  a  bubble  is 

changing  according  to  KXP (- ( (Z-Z01/AL) **2) 

BFT  and  ZOO  -  weighting  function  at  the  bottom  of  a  bubble 

is  changing  according  to  F.XP  (- (  (Z-7.00)  /BET)  **2. 
For  Z.GR.ZO  and  Z.LE.ZOO  the  weighting  function 
is  constant  and  equal  1. 

height  of  the  F-poak  measured 

electron  density  at  the  F-poak  in  e lect rons • cm~ 
electron  density  at  Z*0  in  electrons *cm“ * . 

Note  that  Z»0  corresponds  to  a  height  at  which  compu¬ 
tations  start.  All  distances  are  in  meters. 

All  input  parameters  are  read  from  the  file  DA2.  Hori¬ 
zontal  variations  in  a  medium  or  on  a  phase  screen  are  given 


ZMAX 

F.LMAX 

F.LZ 
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in  arbitrary  units  and  road  from  the  file  ELFL.  Results 
of  computations  are  written  on  the  files  RESULT  and  TERM. 
RESULT  contains  high-accuracy  results  intended  to  be  used 
for  plotting  or  further  calculations.  In  the  later  case 
this  file  must  be  renamed  to  PASS  and  will  be  read  during 
the  next  program  run.  TERM  is  used  only  for  printing  pur¬ 
poses  and  has  the  same  content  as  RESULT. 

Instead  of  using  the  COMPLEX  variables  the  program  uses 
the  matrix  representation  of  complex  variables  with  the 
first  column  being  the  real  part  and  the  second  column  the 
imaginary  part  of  a  variable. 

Subroutine  YO  sets  the  initial  values  to  the  matrix 
representing  the  complex  amplitude.  If  the  phase  screen 
model  is  used,  subroutine  Y01  sets  values  of  the  complex 
amplitude  on  the  screen.  If  KSK1P.NE.0  initial  values  are 
read  f  rom  PASS . 

Using  the  model  of  irregularities,  subroutine  FLUCT 
calculates  dielectric  permittivity  fluctuation  and  its 
mean  value  at  two  adjacent  z  levels. 

Both  A  and  R  matrices  of  equation  (4.7)  are  stored 
in  the  form  of  three  mat  rices:  A,  B,  and  C,  and  AA,  BB  and 
CC.  Matrices  A  and  B  contain  off-diagonal  elements  of  the 
matrix  A,  and  matrix  C  contains  diagonal  elements  of  A. 

Note  that  A~B.  Similarly,  AA  and  BB  contain  off-diagonal 
elements  of  the  matrix  B,  while  in  CC  its  diagonal  elements 
are  stored.  Again  AA=BB .  The  program  was  written  assuming 
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that  A=B  and  AA=BB  but  it  should  not  bo  any  problem  with  its 
general izat ion. 

Elements  of  the  matrix  C  are  calculated  in  the  subroutine 
CTAB  using  output  from  the  subroutine  FLUCT.  If  the  mo<lel 
assumes  no  changes  of  the  electron  density  in  the  z  direc¬ 
tion,  matrices  C  and  CC  can  be  calculated  outside  the  loop 
IX)  2  J-l.N. 

Subroutine  MOZAX  performs  the  matrix  multiplication. 

It  multiplies  the  matrix  A  by  the  matrix  Y  which  contains 
the  complex  amplitude  at  the  preceding  level  or  the  initial 


va 1 ues . 

Subroutine  PROGON  solves 
the  Thomas  algorithm  described 


he  system  of  equations  using 
i n  Sect  ion  4 . 
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High  resolution  electron  density  profile  measured 
inside  thr  plasma  bubble.  (Figure  7  of  McClure  et 
al .  (1977)  reconstructed  after  digi t i ?at ion . ) 


Fig.  1. 


It*)  ItKl 
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Contour  plots  of  constant  AN/N  .it  different  stages 
of  the  bubble  evolution  as  given  by  numerical  sitnu 
lation  (Ossakow  et  a  1 . .  1979). 
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•jq.  4A.  Hori7ont.il  electron  density  fluctuation  profiles  at 
different  heiqhts  inside  a  plasma  bubble,  "initial 
st  aqe . " 
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Fiq.  in.  Horizontal  electron  density  fluctuation  profiles  at 
different  heights  inside  a  plasma  bubble,  "developed 
stage . " 
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Spectral  amplitude  of  small-scale  irregularitie 


Plots  of  the  real  and  imaginary  components  of  the 
complex  amplitude,  amplitude  and  phase  (in  n  units) 
of  a  wave  passinq  throuqh  the  initial  staqe  bubble 
and  received  on  the  qround.  Tin*  wave  frequency  is 
1500  MHz.  Distance  is  measured  from  the  projection 
of  the  center  of  the  bubble  on  the  qround. 


Fiq.  BA.  The  same  plots  as  in  F inure  6  but  .it  136  MHz 

frecpienry.  Th^  wave  structure  at  300  km  height 
is  given. 
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Fiq.  RH .  The  same  plots  as  in  Figure  6  but  at  136  MHz 
frequency.  The  wave  structure  on  the  ground 
is  given. 


OISTancC (KM) 


g .  OB.  W.ivc  structure  after  passing  throuqh  a  phase  screen 
with  discontinuity  at  a  distance  3  km.  The  phase 
screen  height  is  3r>0  km  at  360  MHz. 


Comparison  of  the  amplitude  and  phase  pattern  at 
360  MHz  produced  by  a)  a  phase  screen  and  b)  a 
bubble . 


62 


*» ) 


DISTANCE l«w I 


b) 


Fiq.  11 


The  same'  as  in  Fiquro  10  but  at  4000  MHz  frequency 
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x.*  T 


1. 


4000  MHz 


c.  3 


1 36  MHz 


Distance  Hen) 

Fiq.  1?.  Amplitude  pattern  at  different  frequencies  produr  d 
by  the  "initi.il  staqe"  bubble. 


OIst  »nci  i  ten ) 


Fi<).  1 S .  Wave  structure  at  116  MHz  produced  by  strong 

irreqularit les.  Note  lower  level  of  the  signal 
close  to  the  center  of  the  bubble. 
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Fig.  16A.  Horizontal  electron  density  variations  obtained  by 
randomizing  phases  of  the  Fourier  analyzed  high 
frequency  component  of  the  data  shown  in  Figure  1. 
Fun  A . 
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16n.  Horizontal  electron  density  variations  obtained  by 
randorei z  inq  phases  of  the  Fourier  analyzed  high 
frequency  comf>onont  of  the  data  shown  in  Figure  1. 
Run  H . 
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Amplitude  scintillation  pattern  at  360  Mllz  and 
4000  MHz  produced  by  the  "initial  staqe"  bubble 
with  the  horizontal  structure  shown  in  Fiqurc  16. 
Two  upper  panels  correspond  to  run  A  and  two  lower 
panels  to  run  B. 


